function iq = RFtoIQ2(data,fs,fc,B)
%complex demodulation
%iq=RFtoIQ(data,fs,fc,B)
%data: vetor de dados
% fs : frequencia de amostragem
% fc : frequencia central do transdutor
% B : Banda do transdutor

% retira a parte negativa do espectro 
sp=hilbert(data);
%down mixing
t=(0:1/fs:size(data,1)*(1/fs)-(1/fs))';
t=repmat(t,1,size(data,2));
z=sp.*exp(-1i*2*pi*fc*t);
%low pass
[b,a]=butter(4,fc*B*0.5/fs);%butterworth bandpass filter
s=filter(b,a,z);
%preserva a energia do sinal
s=sqrt(2)*s;
%decimation
dec=4;%
%$dec=ceil(fs/(B*fc));
iq=s(1:dec:end,:);
%iq=s;


